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ABSTRACT 

We studied the X-ray dust scattering halos around 17 bright X-ray point 
sources using Chandra data. All sources were observed with the Advanced CCD 
Imaging Spectrometer and High-Energy Transmission Grating in Continuous 
Clocking Mode (CC-mode) or Timed Exposure Mode (TE-mode). We use an 
iterative method to resolve the halos at small scattering angles from the zeroth 
order data in CC-mode or the hrst order data in TE-mode which is not or less 
piled-up. Using the interstellar grain models of Weingartner & Draine (2001, 
WDOl) and Mathis, Rumpl & Nordsieck (1977, MRN ) to £t the halo prohles, 
we get the hydrogen column densities and the spatial distributions of the scat¬ 
tering dust grains along the line of sights (LOS) to these sources. We hnd that 
the scattering dust density very close to these sources is much higher than the 
normal interstellar medium. For X-ray pulsars GX 301-2 and Vela X-1 with com¬ 
panions of strong stellar winds, the X-ray absorption column densities are much 
higher than the derived scattering column densities, because of the dense media 
around the X-ray sources produce extremely small angle scatterings which are 
indistinguishable from the point sources even with ChandraA angular resolution. 
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For LMC X-1, most of the scattering and absorption occur in Large Magellanic 
Cloud, rather than in the Milky Way. From the obtained X-rays spectra, the 
cold gas absorption and thus the equivalent hydrogen column is determined. We 
have obtained the linear correlation between Nh derived from spectral hts and 
the one derived from the grain models WDOl and MRN (except for GX 301- 
2 and Vela X-1): Nh,wdoi = (0.720 ± 0.009) x NH,abs + (0.051 ± 0.013) and 
Nn,mrn = (1-156 ± 0.016) X Nff^abs + (0.062 ± 0.024) in the units 10^^ cm“^. 
High angular resolution X-ray dust scattering halos offer an excellent potential 
for probing the spatial distributions of interstellar dust medium. 

Subject headings: dust — scattering — X-rays: ISM 


1. Introduction 

X-rays are absorbed and scattered by dust grains when they travel through the inter¬ 
stellar medium. The scattering within small angles results in an X-ray “halo”. The halo 
properties are signihcantly affected by the energy of radiation, the optical depth of the scat¬ 
tering, the grain size distributions and compositions, and the spatial distribution of dust 
along the line of sight (LOS) (Overbeck 1965, Martin 1970, Catura 1983). Therefore an¬ 
alyzing the the X-ray halo properties is an important tool to study the interstellar grains, 
which play a central role in the astrophysical study of the interstellar medium, such as the 
thermodynamics and chemistry of the gas and the dynamics of star formation (Draine, 2003). 

Overbeck (1965) was the hrst to discuss the existence of dust scattering halos. Mathis, 
Rumpl & Nordsieck (1977) derived the dust grain size distributions assuming a silicate- 
graphite dust mixture. Rolf (1983) hrst observationally conhrmed the existence of halos 
using the the data of GX 339-4 observed with the Einstein Observatory. Predehl & Schmitt 
(1995) studied the X-ray halos by systematically examining 25 point sources and 4 supernova 
remnants with ROSAT observations. They found a good correlation between the simulta¬ 
neously measured dust and hydrogen column densities, as well as a correlation between the 
visual extinction and the X-ray derived dust scattering optical depth. Because of the poor 
angular resolution of the ROSAT, the halos in small angles, e.g. less than 50 arcsec, cannot 
be determined accurately; the hat angular distribution of X-ray halos in small angles de¬ 
termined with ROSAT may be due to the underestimation of the halo intensity in the core 
region. Draine & Tan (2003) analyzed the halos of X-ray Nova V1974 Gygni 1992 with the 
ROSAT observations to test the interstellar dust model WDOl and found that the WDOl 
model is very consistent with the observed X-ray halos. With the high angular resolution, 
good energy resolution and broad energy band, the Chandra AGIS is so far the best instru- 
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ment for studying the X-ray halos. Tan & Draine (2004) considered arcsecond scale X-ray 
scattering halos from Galactic Center sources with the Chandra observations and showed 
the halo prohle changed with the spatial distribution of the dust that is close to the source. 
Smith, Edgar & Shafer (2002) presented the observations of the X-ray halo around the Low 
Mass X-ray Binary (LMXB) GX 13-1-1. However, the direct images of bright sources ob¬ 
tained with Advanced CCD Imaging Spectrometer (ACIS)^ usually suffer from severe pileup. 
They therefore investigated only the halo in off-axis angles greater than 50" to avoid any 
possible pileup contamination. 

Making use of the assumption that the real halo could be an isotropic image, Yao et al 
(2003, here after Paper I) reported the reconstruction of the images of X-ray halos in small 
angles from the zeroth order data in CC-mode or the hrst order data in TE-mode. Xiang, 
Zhang & Yao (2004, here after Paper II) improved the method by solving the equations with 
iterative method and using the observational data to generate the Chandra Point Spread 
Function (PSF). They applied the method to the bright sources Cygnus X-1 and Cygnus 
X-3 and fonnd that the halos in small angles are not flat, which is different from the result of 
Predehl & Schmitt (1995). Since the spatial distribution of the dust grains, especially very 
close to the sonrce, is sensitive to the core region of an X-ray halo (Mathis & Lee 1991), it 
is necessary to stndy systematically the halos in small angles for many X-ray point sonrces. 

In this paper, we apply our method described in Paper II, to systematically study the 
X-ray halos of 17 X-ray point sonrces. We resolve the point source halo in very small angles 
from the zeroth order data in CC-mode or the grating data in TE-mode with Chandra, and 
use the different interstellar grain models to fit the halos and get the spatial distribntion of 
the interstellar dust medium (§2). We discuss our results in §3 and present onr conclusions 
^ 4 . 


2. Observations and data reduction 

In this study, we utilized 31 Chandra ACIS High Energy Transmission Grating (HETG)^ 
observations on total 19 sources and all the data are available in the archive npto November 
2004. All these sources (except for PKS 2155-304 and Her X-1) are bright X-rays sonrces 
and have been observed for at least 8.9 ks, therefore there are enongh photons (abont > 0.5 
million connts) for us to resolve the X-ray halos. Both PKS 2155-304 and Her X-1 are very 


^http://cxc.harvard.edu/proposer/POG/html/ACIS.html 
^http://space.mit.edu/HETG/index.html 
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lightly-absorbed X-ray sources with the column density Nh < 3 x 10^'^ cm“^. We use these 
two sources to generate the PSF. We list the source properties and the observation logs in 
Table 1. 


2.1. Halo reconstruction 


Although the zeroth order data in CC-mode and the hrst order data in TE-mode have 
either no or less serious pileup, the zeroth order data in CC-mode have only one dimensional 
images and the hrst order image in TE-mode are mixed with photons of different energies 
and radii, from which we cannot get the halo’s radial prohle directly. Making use of the as¬ 
sumption that the real halo could be an isotropic image, we have reported the reconstruction 
of the images of X-ray halos from the data obtained with the HETGS and/or in CC-mode. 
The method has been described in detail in Paper I and Paper II. Here, we briehy review 
the method and process of the halo reconstruction. 

If the hux of a point source plus its X-ray halo is isotropically distributed and centered 
at the point source as E(r), and the projection process, in which the two dimensional halo 
image is projected to one dimensional image, can be represented by a matrix operator M(r, 
d), then the projected flux distribution P(d) is 

P{d) = M{r,d) X F{r), (1) 


where r is the distance from the centroid source position and d is the distance from the 
projection center (refer to Eig. 1. cited from Paper II). In CC-mode, we can only get the 
count rate C{d), but not the flux projection P{d) directly. With the exposure map of CCDs 
calculated, we can get another equation 

X M(r, d) X F{r) = M'{r, d) x F{r) = C'(d), (2) 


/exposure map 
y matrix 


where M'{r,d) = ^ M{r,d) is another matrix. Using the steepest descent 

method (Marcos & Benar 2001), we can solve equation 1. The iterative process can be 
expressed as 


p{k+l) 




[P-MFW]^[P-MFW] ^ 

[M(P-MPU))]'r[p_MPU)] ^ 


( 3 ) 


where pU+^) and P^^^ are the values of E(r) in the {k + l)th and kih. iterative loops, [P — 
MpWjr 

is the transpose of the matrix P — MF^^k In our iterative process, the loop is 
stopped when A < 0.05 or A reaches its minimal value, where 



-5- 


Table 1. List of sources (column 1) with the observation IDs (column 2),observation mode 
(column 3), positions in equatorial (column 4, 5) and galactic coordinates (column 6 , 7), 
the total counts in 1.0-5.0 keV (column 8 ), and the total exposure time (column 9) 


Source 

ObsID 

mode 

R. A. [2000] 

Dec. [2000] 



counts 

time (ks) 

4U1705-44 

1923 

TE 

17'*08’"55'’ 

-44‘*06'oo" 

343.3"^ 

-2.3'* 

2637040 

25.3 

4U1728-16 

703 

TE 

17'‘31™44'’ 

-16‘*57'42" 

00 

cn 

9.0'* 

2619501 

20.9 

Cir X-1 

1700, 1905 
1906, 1907 

TE 

15/i20™41'’ 

-57‘*10'0l" 

322.1'* 

0 

d 

5208035 

40.3 

Cyg X-1 

107, 3814 

TE 

19'‘58’"22s 

+35'*12'o6" 

71.3'* 

3.1'* 

8225397 

58.6 

Cyg X-2 

1102 

TE 

2l'i44™4is 

+38^*19'is" 

87.3'* 

-11.3'* 

3513440 

14.6 

Cyg X-3 

425 

TE 

20'*32’"27s 

+40‘*57'l0" 

79.8d 

0.7'* 

1431689 

16.2 

GRS 1915-1-105 

660 

TE 

19'*15’"12s 

+10‘*56'44" 

45.4'* 

0.2'* 

3379930 

30.6 

GX 13-1-1 

2708 

TE 

18'*14"*32s 

-17^*09' 27 " 

13.5'* 

0.1'* 

2857022 

39.9 

GX 301-2 

2733, 3433 

TE 

12'*26"*38s 

-62‘*46'l3" 

300.1'* 

0.0'* 

584812 

98.2 

GX 3-1-1 

2745 

GC 

17'‘47™56'’ 

-26^*33' 49 " 

2.3'* 

00 

d 

1298798 

9.2 

GX 340+0 

1921 

TE 

16'*45"*48s 

-45^*36'40" 

339.6'* 

-0.1'* 

1385351 

24.0 

GX 349+2 

3354, 715 

TE 

17'‘05™45® 

-36^*25' 23 " 

349.1'* 

2.7'* 

6602027 

27.0 

GX 5-1 

716 

TE 

18'‘01’"08'’ 

-25^*04' 45 " 

5.1'* 

-1.0'* 

2123066 

89.1 

GX 9+1 

717 

TE 

18'‘01™32'’ 

-20'*3l'44" 

9.1'* 

1.2'* 

1275273 

89.9 

LMG X-1 

93 

TE 

05/139171408 

-69^*44' 37 " 

280.2'* 

-31.5'* 

455010 

21.8 

Ser X-1 

700 

TE 

18'‘39’"58'’ 

+05'*02'09" 

36.1'* 

4.8'* 

7171223 

76.3 

Vela X-1 

1927, 1928 

TE 

09/io2™07'’ 

-40‘*33'l7" 

263.1'* 

3.9'* 

1374192 

62.1 

Her X-1 

2749 

TE 

16/157771508 

35 ^* 20 ' 33 " 

58.1'* 

37.5'* 



PKS 2155-304 

1014,1705,3167 

TE 

2l'‘58’"52'’ 

-30‘*13'32" 

17.7'* 

-52.2'* 




333 , 3706,3708 
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N is number of P{d) and AP is the error of P{d). Replacing M in equation 3 with M', as 
well as P with C, we can solve equation 2. 

The accurate Chandra PSF is important for reconstructing the halo accurately. Follow¬ 
ing the same as procedure described in Paper II, we calculate the zeroth order PSF from 
the observation data of bright point sources without halo (for example Her X-1 and PKS 
2155-304) at large angles and from the MARX^ simulation at small angles. Because the data 
of Her X-1 and PKS 2155-304 in HETGs are not suffered from pileup even when the source 
is pointed on-axis, the projection of grating arms PSF in small angles is also calculated from 
the observation data. In order to improve the statistical quality, we use data from seven 
observations to generate the PSF. These data IDs are listed in Table 1. 

For the TE-mode data, we use the hrst order data (HEG ±1 and MEG ±1) within 
60 arcsec around the source position. First, we divided the dataset into 20 energy bands, 
spaced every 0.2 keV from 1.0-5.0 keV, and use the process described in Paper H to generate 
the pure halo projections. Then we sum them from each energy bands and derive the halo 
radial prohle in units of photons cm“^ s“^ arcsec“^. For zeroth order of the GG-mode data, 
we also use the same process of data reduction as described in Paper H to obtain the pure 
halo radial flux in each energy band. Then we sum the halo radial flux from each energy 
band and get the total halo radial flux from 1.0-5.0 keV. 


2.2. Halo model fittings 

As discussed by Mathis & Lee (1991), the geometry of the X-ray scattering can be 
demonstrated as Figure 2. Because all scattering angles are small, 6* ~ (1 — x)6sca is always 
a good approximation, where x = d/D, d is the distance of the scattering grain from the 
observer, and D is the distance to the source. 

The observed halo surface brightness lsca{0) at an angle 6 from the point source of 
X-rays can be shown as 

^ = Nh dES{E) dan{a) f dxf{x){l - x)-^ x ^ ( 4 ) 

Jeo Ja^i„ Jo dil 

where Nh is the total column density of hydrogen between the observer and the X-rays 
source, n{a)da is the number of grains per H atom with radii between a and a + da, Ex is 
the observed flux, S{E)dE is the normalized photon energy distribution of the X-ray point 


^http://space.mit.edu/ASC/MARX/ 
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Zeroth Order data 



EO 


P(d) 



Fig. 1.— The projection of the photons in zeroth order image along the read-out direction 
and the projection of the photons along a grating arm. 



Fig. 2.— Geometry of the X-ray scattering process. 




source, f{x) is the density of hydrogen at xD relative to the average density along the line 
of sight to the X-ray source, and is the differential scattering cross section for a single 
grain. 

Since the X-ray halo is determined by the size, position and composition of the dust 
grains, as well as the source flux and absorption column, we need to make some assumptions 
about the dust grains before we can use equation 4 to £t the halo profiles. We used the 
grain models Weingartner & Draine (2001, WDOl) and Mathis, Rumpl & Nordsieck (1977, 
MRN) to fit our halo radial prohle of these sources. The MRN model is a classical model 
with the size distributions of dust grains: n(a) oc (0.005 fj,m < a < 0.25 jj,m). The size 

distribution of dust grains in WDOl model is rather complicated, which includes sufficient 
very small carbonaceous grains and larger grains (please refer to Weingartner & Draine 2001 
for detail). Since we only concern the halo in small angles (< 60”, except for the source 
LMC X-1), multiple scatterings can be neglected even though the scattering optical depth 
Tsca is as large as 2. The model htting codes were provided by Randall K. Smith and the 
same codes have been applied to the halos of GX 13-1-1 (Smith, Edgar & Shafer 2002). 

As a first try, we use the smoothly distributed dust models to fit the halo. We find that 
these distributions cannot describe the whole halo profile of any source we analyzed, - all the 
xi are greater than 5.0 for both the WDOl and MRN models. Then we used the “extended” 
MRN model (Witt, Smith & Dwek 2001) which has the same total mass as the MRN model, 
but extends the MRN size distributions to 2.0 fim to fit the halo radial profile. We find 
that these models underestimate the photons in large angles (>20 arcsec) and in very small 
angles (about 1 3 arcsec), and overestimate the photons in 3 20 arcsec. The same is true for 
the ’’extended” WDOl model which extends to the size distributions to 2.0 to 5.0 /rm. We 
then divide the LOS to four parts (a;:0.0-0.25, 0.25-0.5, 0.5-0.75, 0.75-1.0), in order to probe 
the inhomogeneity at large scales, and assume that the dust grains are smoothly distributed 
in each part and that the amount of dust in each part is allowed to vary independently. 
However we found that there are systematic excess halo intensity in small angles (below 10 
arcsec), indicating that there should be substantial dust grains at the position very close 
to the X-ray sources. We also found that the fitting are not very sensitive in constraining 
the relative amount of dust grains located at 0.0-0.25 and 0.25-0.50, but are very sensitive 
in the one located at 0.75-1.0. We therefore re-divide the LOS to logarithm-equally spaced 
bins with more bins close to the sources: 0.0-0.651, 0.651-0.881, 0.881-0.962, 0.962-0.99 and 
0.99-1.0; in this scheme the amount of dust grains in each part may be determined relatively 
independently. Last we try to use the power-law or exponential function to fit the dust 
spatial distributions and find the fit results are very poor. The best fit results and the 
values of for model are listed in Table 2 and Table 3 for the WDOl and MRN models, 
respectively. The X-ray halo radial profiles and relative density distributions along the LOS 
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for all sources are shown in Figure 3. The fractional halo intensity is calculated from each 
halo from l.S" to 60". 

The source LMC X-1 is treated differently from other sources, since it is not in our 
galaxy. The detailed htting process is described in section 3. The halo of LMC X-1 from 
l.O" to 40” is obtained from the grating data and that from 40" to 180" is extracted directly 
from the zeroth order ACTS image. 


2.3. Spectral Fittings 

Here we perform spectral hts to derive the hydrogen absorption column densities with 
the high resolution spectra from these sources observed with ACIS-S -|- HETG, in order to 
compare with the scattering column densities derived from X-ray halos. We use the ciao 
3.0.2 to extract the grating spectra. In order to reduce the contamination of the halo, the 
region to hlter the data is limited in 5 arcsec for the TE-mode data. 

The HEG±1 and MEG±1 spectra sometimes suffered from pileup when the energy is 
higher than 2 keV because some sources are very bright, i.e, Cygnus X-1 and Circinus X- 
1. However the energy spectra from MEG±3 are not affected by the pilenp and there are 
enongh photons for spectral analysis. Althongh the third order photons below 2 A may 
be contaminated by pilenp from the MEG hrst order for the brightest source, e.g, Circinus 
X-1 (Schulz & Brandt, 2002), this has no impact on our analysis because we only use the 
spectra below 5 keV (~ 2.4 A). Since there are less photons in the low energy bands (e.g. 
less than 2 keV) of MEG±3 spectra and sometimes the MEG±3 spectra are cnt off below 2 
keV, the HEG±1 or MEG±1 energy spectra below 2 keV band and the MEG±3 spectra in 
2-5 keV bands are also used. The spectra are rebinned so that there are at least 100 connts 
in each energy bin for the MEG±3 spectra, as well 400 connts in each energy bin for the 
MEG±1 or HEG±1 spectra. We use the Xspec 11.1.0 to ht these energy spectra. Three 
models, i.e., the power law, black body and thermal bremsstrahlnng, are used, including cold 
gas absorption, which was modeled using the absorption cross sections given by Morrison & 
McCammon (1983). For each source, we only list the results for one model which can best 
fit the data. The highly accurate emission lines in these spectra are ignored since we only 
concern the absorbing hydrogen column density. A snmmary of our spectral htting results 
are given in Table 4. 
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Fig. 3.— The X-ray halo prohles of the point X-ray sources, htted using the interstellar 
grain models WDOl and MRN. The solid lines are the model WDOl and the dashed lines 
are the model MRN. The relative hydrogen density spatial distributions are shown as insets 
in the figure, where the horizontal axis is for the relative distance x to the source and the 
vertical axis is the relative density normalized to the value in the last bin (close to 1). 

Table 2. The relative hydrogen densities in different segments along the LOS (column 2-6, 
relative to the peak hydrogen density along each LOS) and total hydrogen column densities 

(column 7) from the grain model WDOl htting. 


Source 

The relative hydrogen density in different position from WDOl fit 
0.0-0.651 0.651-0.881 0.881-0.962 0.962-0.99 

0.99-1.0 

total Nh 
(10^^ cm“^) 

X^/dof 

4U1705-44 

0.0331 ±0.0019 

0.0118 ±0.0020 

0.0000 ± 0.0034 

0.027 ±0.010 

1.00 ±0.08 

0.964 ± 0.044 

13.1/24 

4U1728-16 

0.0201 ±0.0031 

0.0000 ± 0.0029 

0.019 ±0.005 

0.019 ± 0.013 

1.00 ±0.08 

0.241 ±0.022 

4.1/24 

Cir X-1 

0.0672 ± 0.0019 

0.0138 ± 0.0020 

0.0359 ± 0.0035 

0.013 ± 0.009 

1.00 ±0.08 

1.77 ±0.05 

32.3/24 

Cyg X-1 

0.000 ± 0.005 

0.0807 ± 0.0034 

0.106 ±0.005 

0.092 ± 0.013 

1.00 ±0.07 

0.413 ±0.021 

41.0/24 

Cyg X-2 

0.0063 ± 0.0019 

0.0033 ±0.0017 

0.0000 ± 0.0031 

0.029 ± 0.008 

1.00 ±0.05 

0.295 ±0.027 

4.1/24 

Cyg X-3 

0.0587 ± 0.0020 

0.0167 ± 0.0024 

0.0115 ±0.0043 

0.051 ± 0.012 

1.00 ±0.13 

2.49 ± 0.09 

17.4/24 

GRS 1915±105 

0.217 ±0.005 

0.119 ±0.005 

0.036 ± 0.009 

0.244 ± 0.023 

1.00 ±0.21 

2.85 ±0.06 

33.3/24 

GX 13±1 

0.0590 ± 0.0025 

0.0916 ±0.0029 

0.065 ± 0.005 

0.091 ± 0.013 

1.00 ±0.11 

2.08 ±0.06 

44.2/24 

GX 301-2 

0.0091 ± 0.0009 

0.0013 ±0.0012 

0.0006 ± 0.0022 

0.004 ± 0.008 

1.00 ±0.13 

3.55 ±0.33 

3.8/24 

GX 3±1 

0.107 ±0.005 

0.108 ±0.006 

0.212 ±0.015 

0.48 ± 0.07 

1.00 ± 0.74 

1.35 ±0.09 

40.9/24 

GX 340±0 

0.0710 ±0.0012 

0.0029 ± 0.0015 

0.0250 ± 0.0026 

0.044 ± 0.008 

1.00 ±0.14 

4.82 ± 0.08 

19.2/24 

GX 349±2 

0.0310 ±0.0023 

0.0127 ± 0.0022 

0.0000 ± 0.0041 

0.068 ± 0.012 

1.00 ±0.08 

0.611 ±0.032 

8.8/24 

GX 5-1 

0.0439 ± 0.0018 

0.0242 ± 0.0021 

0.0532 ± 0.0040 

0.128 ± 0.011 

1.00 ±0.10 

2.00 ± 0.06 

10.5/24 

GX 9±1 

0.0948 ± 0.0038 

0.0041 ± 0.0040 

0.000 ± 0.008 

0.138 ± 0.020 

1.00 ±0.16 

1.074 ± 0.044 

9.9/24 

Ser X-1 

0.0278 ± 0.0028 

0.0000 ± 0.0027 

0.019 ±0.005 

0.015 ± 0.012 

1.00 ±0.07 

0.329 ± 0.023 

5.2/24 

Vela X-1 

0.0158 ±0.0031 

0.0186 ± 0.0032 

0.0037 ± 0.006 

0.034 ± 0.016 

1.00 ±0.11 

0.436 ± 0.041 

3.2/24 


x(0.0 - 1.0) : 0.0 - 0.737^9 - 0.995 0.9957^019 “ 


LMC X-1 relative density 0.0142 ± 0.0017 0.0033 ± 0.0003 1.00 ±0.08 




0.444 ± 0.027 35.3/41 
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Table 3. The lists of relative hydrogen densities in different positions (column 2-6, relative 
to the max hydrogen densities in this sight) and total hydrogen column densities (column 

7) from the grain model MRN £t. 


Source 

0.0-0.651 

Nh (10^^ cm“^ 
0.651-0.881 

) in different position from MRN fit 
0.881-0.962 0.962-0.99 

0.99-1.0 

total Nh 
(10^^ cm~^) 

X^/dof 

4U1705-44 

0.0382 + 0.0031 

0.0276 + 0.0026 

0.0052 + 0.0041 

0.000 + 0.010 

1.00 

+ 

0.07 

1.51 + 0.08 

13.0/24 

4U1728-16 

0.023 ± 0.006 

0.0098 + 0.0041 

0.015 + 0.006 

0.014 + 0.012 

1.00 

+ 

0.08 

0.362 + 0.047 

5.7/24 

Cir X-1 

0.0838 ± 0.0032 

0.0408 + 0.0028 

0.0495 + 0.0043 

0.005 + 0.010 

1.00 

+ 

0.07 

2.91 + 0.09 

33.2/24 

Cyg X-1 

0.000 ± 0.007 

0.0546 + 0.0041 

0.165 + 0.006 

0.056 + 0.012 

1.00 

+ 

0.06 

0.596 + 0.048 

34.8/24 

Cyg X-2 

0.0037 ± 0.0031 

0.0094 + 0.0022 

0.0000 + 0.0030 

0.000 + 0.007 

1.00 

+ 

0.05 

0.38 + 0.06 

5.2/24 

Cyg X-3 

0.0662 ± 0.0031 

0.0408 + 0.0030 

0.021 + 0.005 

0.023 + 0.013 

1.00 

+ 

0.10 

4.03 + 0.15 

14.4/24 

GRS 1915+105 

0.192 + 0.006 

0.207 + 0.006 

0.054 + 0.009 

0.206 + 0.021 

1.00 

+ 

0.15 

4.59 + 0.11 

32.3/24 

GX 13+1 

0.0581 + 0.0037 

0.1103 + 0.0035 

0.114 + 0.006 

0.070 + 0.013 

1.00 

+ 

0.08 

3.25 + 0.11 

44.0/24 

GX 301-2 

0.0138 + 0.0018 

0.0067 + 0.0020 

0.0000 + 0.0028 

0.000 + 0.008 

1.00 

+ 

0.12 

4.47 + 0.37 

3.6/24 

GX 3+1 

0.070 + 0.006 

0.136 + 0.006 

0.160 + 0.012 

0.468 + 0.049 

1.00 

+ 

0.35 

2.04 + 0.10 

41.2/24 

GX 340+0 

0.1015 + 0.0025 

0.0511 + 0.0024 

0.0147 + 0.0040 

0.093 + 0.010 

1.00 

+ 

0.08 

7.56 + 0.16 

15.0/24 

GX 349+2 

0.0307 + 0.0035 

0.0266 + 0.0027 

0.0005 + 0.0042 

0.028 + 0.011 

1.00 

+ 

0.06 

0.94 + 0.06 

9.3/24 

GX 5-1 

0.0460 + 0.0027 

0.0408 + 0.0026 

0.0595 + 0.0045 

0.127 + 0.011 

1.00 

+ 

0.07 

3.12 + 0.11 

9.6/24 

GX 9+1 

0.097 + 0.005 

0.0353 + 0.0045 

0.000 + 0.008 

0.080 + 0.018 

1.00 

+ 

0.11 

1.83 + 0.08 

9.8/24 

Ser X-1 

0.044 + 0.005 

0.0000 + 0.0035 

0.025 + 0.005 

0.000 + 0.015 

1.00 

+ 

0.07 

0.586 + 0.049 

7.8/24 

Vela X-1 

0.014 + 0.005 

0.0285 + 0.0042 

0.013 + 0.007 

0.003 + 0.017 

1.00 

+ 

0.09 

0.62 + 0.08 

5.5/24 


x(0.0 - 1.0) : 

0.0 - 0.277o;27 

0.807o;o7 - 0.9973 

0.99737“;^^;^ - 1-0 






LMC X-1 

relative density 

0.0043 + 0.0006 

0.0027 + 0.0003 

1.00 + 0.08 




0.75 + 0.05 

34.8/41 


Table 4. 


List of the spectral fits (all Nh values are in units 10^^ cm ^). 


Our fits used 


power-law (pwl), thermal bremsstrahlung (thb) or black body (bb) model. 


Source 

Nh (Lit) 

Nh 

r (pwl) 

kT (thb) 

kT (bb) 

XS/N 

4U1705-44 

1.23L 1.45^, 0.85® 

1.23 + 0.04 



1.01 + 0.02 

0.86/1551 

4U1728-16 

0.26L 0.26^, 0.21® 

0.24 + 0.02 

-1.76 + 0.04 



0.95/1547 

Cir X-1 

2.4lh 1.98® 

2.41 + 0.05 


8.52R;®4 


0.92/574 

Cyg X-1 

0.41L 0.62®, 0.72® 

0.441 + 0.017 

-1.40 + 0.03 



1.61/2029 

Cyg X-2 

0.25L 0.28^, 0.22® 

0.217 + 0.035 

-1.43 + 0.05 



1.27/382 

Cyg X-3 

3.3lh 1.42® 

3.20 + 0.10 



1.51 + 0.06 

1.38/713 

GRS 1915+105 

3.5^, 1.76® 

3.852 + 0.034 



0.900 + 0.010 

0.86/3169 

GX 13+1 

2.79h 1.79® 

2.671 + 0.027 



-L-nu_Q 

0.92/2654 

GX 301-2 

85®, 1.81® 

16.2 + 0.8 

0.13 + 0.12 



0.82/605 

GX 3+1 

1.421, 2® 

1.56 + 0.06 



0.72 + 0.05 

1.05/926 

GX 340+0 

5.01', 2.188 

6.65 + 0.14 


13.0t®d 


1.10/380 

GX 349+2 

0.961, Q ggS 

0.929 + 0.013 

1.44 + 0.03 



0.96/6080 

GX 5-1 

2.781,3.0®, 0.94® 

2.87 + 0.05 



1.020 + 0.020 

0.85/660 

GX 9+1 

1.871, I. 43 I, 0.91® 

1.339 + 0.040 


1.00 + 0.03 


0.80/456 

LMC X-1 

O. 94 I, 0.07® 

0.653 + 0.013 


2.93 + 0.05 


1.41/801 

Ser X-1 

O. 50 I, 0.47® 

0.487 + 0.017 


c 7 c:0.22 

0 . / U_Q 


1.27/2744 

Vela X-1 

1 - 18®, 0.51® 

1.05 + 0.05 

0.98 + 0.08 



0.99/412 


^The literature values (column 2) are taken from Predehl Sz Schmitt 1995 (1), Schulz 1999 (2), Schulz et al 2002b 
(3), Chapuis Corbel 2004 (4), Watanabe et al 2003 (5), Kretschmar et al 1997 (6), Gilfanov et al 2003 (7) and the 
galactic H I column density (8) from heasarc website “http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pr’. 
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3. Discussion of results 

The main quantities derived from our X-ray analysis for each of our sample sources 
are the equivalent hydrogen absorption column density and the hydrogen scattering column 
density Nh, the spatial distribution of the dust medium and the fractional halo intensity 
IfracO--^" < 0 < 60"). In the following we will discuss these results and investigate the 
correlations between these measurements. 


3.1. Spatial distribution of dust medium 

In section 2, we have shown that we can use our method to resolve the halos as close 
as 1.5 arcsec to the surrounded point sources. These small angle halos could provide tight 
constraints on the different interstellar grain models and on the dust spatial distribution 
especially when the dust grains are very close to the sources. 

Both the WDOl model and MRN model hts indicate that there should be substantial 
amount of dust grains near the sources; these high density dust may be the molecular clouds 
in which the binary systems were formed initially (see section 3.5 for more discussions). 

Next we move on to discuss the results of some sources which are either different from 
the majority of other sources, or our results are not exactly the same as the existing results 
in the literature. 


3.1.1. LMC X-1 

LMC X-1 is located in the Large Magellanic Cloud (LMC), about 55 kpc away from us 
and its galactic coordinates are / = 280.2*^, b = —31.5*^. Predehl & Schmitt (1995) observed 
it with ROSAT, and found fhaio = 2.9% at 1.26 keV, less than our result, fhaio = 4.2% from 
1.5 arcsec to 180 arcsec at 1.79 keV, if fhaio is proportional to E~'^. Predehl & Schmitt 
(1995) estimated that the halo and PSF have equal intensities at about 200" for the LMC 
X-1 observation. But according to our analysis, the halo intensity is concentrated mostly 
in small angles: fhaio{l-0'' < 0 < 60") = 2.9% and //ia«o(60" < 9 < 180") = 1.3%. Thus 
it appears possible that the much coarse angular resolution of the ROSAT has resulted in 
under-estimation of total halo intensity. 

Since the LMC X-1 is at high galactic latitude and the X-ray photon travels through 
the two galaxies, we divided the LOS to three parts: the hrst one is in our galaxy (0 ~ xi), 
the second one is in the LMC {x 2 ~ 1) and the last one is nearby the source ( 0:3 ~ 1), where 
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X 2 i X 3 are the positions relative to the distance from LMC X-1. We assume the dust 
grains in each part are smoothly distributed along the LOS. First we change the values of 
Xi, X 2 , X 3 manually to £t the halo profile with the grain model WDOl and get the value of 
for each trial. Searching for the minimal we get the best fit values of Xi^ 2 , 3 - Then we 
separately change the value of xi, X 2 , X 3 and get their uncertainties. The best values are: 
xi = 0.14]'^q;^4, X 2 = 0.73tao9; ^3 = 0-995taooi9- Finally, from the model hts, we get the Nh 
(1022 cm-2); = o'.112, Km = 0-282. The total hydrogen 

column density is 0.444 x IO 22 cm“2 and is consistent with Nh = 0.46 x IO 22 cm“2 derived 
from spectral fitting by Cui et al (2002). For the MRN model, the best-fit parameters are 
similar and listed in table 3. 

These results indicate that the dust grains in the high galactic latitude are much less 
than the ones nearby the LMC X-1 (within about 0.3 kpc of the source), and the dust grain 
density along the LOS (excluding the vicinity of LMC X-1) to the LMC is very low. It is 
consistent with the fact that most dust grains are located in the galactic disk. In fact the 
Nh of LMC X-1 inside the Milky Way derived in our model fittings is in good agreement 
with Nh = 0.07 x IO 22 cm“2^ inferred from the 21-cm survey as listed in Table 4 (indicated 
by literature (8)), which is the HI column density inside the Milky Way but averaged along 
the approximate direction of LMC X-1 (the closest direction with available data is about 0.3 
degree apart from LMC X-1); the small difference may be caused by the inhomogeneity of the 
molecular clouds along the approximate direction of LMC X-1. We mention in passing that 
it is not meaningful to make comparisons between Nh indicated by literature (8) in Table 
4 with Nh derived from either the grain model fittings or spectral httings for all galactic 
sources, because the former is obtained by integrating the HI density through the whole 
Galaxy along the approximate direction of the particular source and is thus more adequate 
for estimating Nh of extra-galactic sources. 


3.1.2. Cygnus X-1 and Cygnus X-3 

Cygnus X-1 is the hrst black hole candidate, located at I = 71.3^, b = 3.1^. It is one of 
the brightest X-ray sources and the estimated distance is D = 2.5 ±0.3 kpc (Ninkov, Walker 
& Yang, 1987). Schulz et al (2002) used the X-ray absorption spectroscopy observed with 
Chandra ACIS-S ± HETG to derive the hydrogen column density, Nh = 6.2 x 102^ cm“2. 
We have derived the fhaio{^-^" < d < 6 O”) = 6.1% at 1.99 keV; this is slightly more than the 
value fhaio = 11.6% at 1.20 keV, found by Predehl & Schmitt (1995), but consistent with 
the result in Paper I and Paper H. 

The grain models WDOl and MRN both show that the dust grains between the observer 
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and Cygnus X-1 are not smoothly distribnted along the LOS and the dnst grains very close 
to the sonrces. This is consistent with the resnlt in the Paper II. However the details of the 
dust spatial distribution derived here is not the same as the result in Paper II, because in 
Paper II we assumed only one segment of scattering dust grains located along the LOS. 

Cygnus X-3 is an X-ray binary with more than 40% halo flux in 0.1-2.4 keV {ROSAT 
energy range) (Predehl & Schmitt, 1995), located at I = 79.8'^, b = O.T^. The bright X-ray 
halo has been used to determine the distance to the source as about 9 kpc (Predehl, et al, 
2000). The halo fraction fhaioi^-^” < 0 < 60”) = 8.5% in 1.0-5.0 keV is much less than the 
total halo fraction of ROSAT which includes the large angle halo, indicating a substantial 
amount of scattering dust is located at the near side to the observer. 


3.1.3. GX 13+1 

GX 13-1-1, a LMXB, is a bright and highly absorbed X-ray source and located at I = 
15.5'^, b = O.l'^. Smith, Edgar & Shafer (2002) observed the GX 13-1-1 with Chandra AClS-l, 
and extracted the halo between 50" and 600”. They found the halo fraction is fhaioi^O" < 
6 < 600”) ~ 12% in 2.8 keV and the total halo intensity as a function of energy is I{E) = 
l.hlgjE^gy and I{E,50” <6 < 600”) = (0.939 ± 0.028)77^^^. This corresponds to fhaio{0" < 
6 < 50”) = (7.02 ~ 13.54)% at 2.8 keV. This is consistent with our result fhaioi^A” < 6 < 
60") = 11.4% at 2.82 keV. 

Our fits show that the total hydrogen column density Nh = 2.08 x 10^^ cm“^ from 
the grain model WDOl is signihcantly larger than the value 1.74 x 10^^ cm“^ derived from 
the WDOl model £t by Smith, Edgar & Shafer (2002). It is possible that Smith, Edgar & 
Shafer (2002) underestimated the dust located at a; > 0.90 which does not affect the observed 
I{9h > 50"), since these dust near the source contributes to the halo mainly in angles smaller 
than 50". This interpretation is also consistent with the hydrogen scattering column density 
Nh = 1.64 ± 0.07 x 10^^ cm“^ derived from the WDOl htting by integrating the hydrogen 
density in a; = 0.0 — 0.90 from the results in Table 2. 


3 . 1 . 4 . GX 301-2 and Vela X-1 

The X-ray binary pulsar GX 301-2 is located at / = 300. l'^, b = 0.0“^. The system 
consists of an accreting magnetized neutron star in a highly eccentric orbit, embedded in 
the stellar wind from a B2 supergiant companion star. When the orbital phase is in the 
pre-periastron, the hydrogen column density can be as large as 1 x 10^^ cm“^ (Watanabe et 
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al 2003). 

The hydrogen column density from the grain model £t is much smaller than the one 
from the spectral fit. A likely interpretation is that there are strong and dense stellar winds 
around the X-ray source and they absorb and scatters the X-ray photons severely. Therefore 
the Nh from the spectral hts reflects the total absorption column density correctly. However 
the scattered photons in the halo caused by the stellar winds are very close to the point 
source as x is very close to 1, e.g., 9 < O.S" for < 0.1 pc which corresponds to 

X > 0.9999, where r^jind is the radius of the dense stellar winds. Consequently this part of 
halo is indistinguishable from the point source, resulting in underestimation of the scattering 
column density. 

Vela X-1, similar to GX301-2, is an eclipsing high-mass X-ray binary consisting of an 
early type supergiant and a pulsar, having an orbital radius of about 1.7/?* (~ 53Rq) with a 
period of 8.964 days. The X-ray pulsar is accreting from the intense stellar winds emanating 
from the companion and the X-ray spectrum shows considerable photoelectric absorption at 
lower energies, varying with the orbital phase, as the neutron star is deeply embedded in 
the dense stellar wind (Haberl & White, 1990). Kretschmar et al (1997) found the column 
density Nh varies from (1 ± 1) x 10^^ cm“^ to (18 ± 4) x 10^^ cm“^. Since the halo caused 
by the stellar winds has an angular size smaller than 0.5", the scattering column density is 
also underestimated. 

For several other sources, e.g., 4U 1705-44, Cir X-1 and Ser X-1, the halo intensity 
between 1.5" and 60" is also smaller than the total halo intensity obtained by Predehl & 
Schmitt (1995), most likely due to the smaller angular size studied here. 


3.2. Correlations between different Nh column densities 

In Figure 4 we show the correlation between Nh,mrn and Nh,wdoi which is well fitted 
by a linear relation, Nh,mrn = (1.594±0.023) x Vj:/vrDoi + (~0.009±0.025). This shows that 
the column density derived by fitting the X-ray halo using the MRN model is ~ 60% larger 
than the fitting results using the WDOl model. This is close to the result (~ 45%) of Smith, 
Edgar & Shafer (2002). The systematic difference for the inferred Nh in the two models 
is most likely due to more large-grains in WDOl model than in MRN model, which scatter 
photons preferentially in small angles, i.e, more total scattering cross section a'^dnia) 

^ ^min 

of the dust per H nucleus in the WDOl model, about 65% more than in the MRN model. 

In Figure 5 we show the correlations between the dust scattering column densities 
Nh,wdoi and Nh,mrn derived from the two models, and X-ray absorbing column density 
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NH,abs derived from X-ray spectral fitting. It is obvious from Figure 5 that both correlations 
are extremely tight, which can be described by two linear relations Nh^wdoi = (0.720 ± 
0.009) X NH,abs + (0.051 ± 0.013) and Nm,mrn = (1-156 ± 0.016) x NH,abs + (0.062 ± 0.024), 
in units of 10^^ cm“^. The outliers, i.e. GX 301-2 and Vela X-1 as discussed earlier, are not 
included in the fitting. 


3.3. FHI and Nh 


The fractional halo intensity (FHI) I/rac can be shown as (Mathis & Lee 1991): 

I halo 


Ifrac 


i^^halo T j 


= (l-e- 


( 5 ) 


where the optical depth in scattering r^ca equals to Tsca{E = IkeV) x {Emean/^keVY within 
Gaussian approximation. The value for the mean energy Emean has been obtained by weight¬ 
ing each X-ray photon energy according to the energy dependence of the scattering cross 
section, i.e. approximately E~‘^. The fractional halo intensity for all sources and the mean 
energy Emean are shown in Table 5. 


Since Tsca{E = IkeV) ~ Nh x {E/lkeV) equation 5 can be rewritten to 

k X N [10^^ cm~^] 

7/mc = l-e ^ , (6) 


where k is approximately a constant. Using the result of Predehl & Schmitt (1995) to 
estimate the value of /c, we can get k ~ 0.5. 

In our analysis we calculated only the halos between 1.5” and 60”. The fractional halo 
intensity //rac(l-5” < 9 < 60") will be less than Ifradtotal). For each source, Ifracd-5" < 
9 < 60”) = m X Ifrac{iotal), m depends on the mean energy, the dust spatial distribution 
and the total hydrogen column density of each source, i.e., m varies across the sources. 
However in Figure 6, a good correlation is quite obvious, indicating that m does not change 
signihcantly for the sources we have analyzed, excluding the two outliers GX 301-2 and Vela 
X-1. The solid curve in Figure 6, which best describes the correlation, is a plot of equation 
6 with m = 0.53. 


3.4. Correlation between Nh and the galactic latitude 

If the scale height of dust grains in the disk of Galaxy is d and the X-ray point source 
is outside the dust grain layer, the effective distance which X-ray photons travel through 
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Nh(WDOI) (1022 cm-2) 

Fig. 4.— Correlation between the scattering column density Nh derived by fitting the X-ray 
halos with MRN and WDOl models respectively. 



Fig. 5.— The correlations between the dust scattering column densities Nh^wdoi and 
Nh,mrn derived from the two models, and X-ray absorbing column density NH,abs derived 
from X-ray spectral htting. 
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Table 5. Lists of the fractional halo intensity (column 2), the mean energy weighted by 
the scattering cross section (column 3) and the hydrogen absorption column densities 

(column 4). 


Source 

Ifrac{%} 

(EmeanllkeV)^ 

Nh(10^^ cm 

Nh 

4U 1705-44 

4.6 

2.53 

1.23 + 0.04 

0.192 

4U 1728-16 

2.1 

1.90 

0.238 + 0.015 

0.066 

Cir X-1 

8.0 

2.64 

2.414 + 0.048 

0.346 

Cyg X-1 

6.1 

1.99 

0.441 + 0.017 

0.111 

Cyg X-2 

2.4 

1.94 

0.217 + 0.035 

0.058 

Cyg X-3 

8.7 

3.20 

3.20 + 0.10 

0.313 

GRS 1915-1-105 

12.3 

2.94 

3.852 + 0.034 

0.446 

GX 13+1 

11.5 

2.82 

2.671 + 0.027 

0.336 

GX 301-2 

5.8 

4.00 

16.2 + 0.8 

1.012 

GX 3+1 

6.9 

2.46 

1.56 + 0.06 

0.258 

GX 340+0 

14.9 

3.25 

6.65 + 0.14 

0.630 

GX 349+2 

3.9 

2.24 

0.929 + 0.013 

0.185 

GX 5-1 

9.7 

2.90 

2.87 + 0.05 

0.342 

GX 9+1 

4.9 

2.55 

1.339 + 0.040 

0.206 

LMG X-1 

2.9 

1.79 

0.653 + 0.013 

0.204 

Ser X-1 

2.5 

1.95 

0.487 + 0.017 

0.128 

Vela X-1 

2.9 

2.35 

1.05 + 0.05 

0.190 



NH(l022cm 2)/E2(keV) 

Fig. 6.— The correlation between the fractional halo intensity (FHI) and the hydrogen 
column density. The solid line is plotted according to equation 6 with m = 0.53. 
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the galactic dust layer is d/sin{\tf^\), where is the galactic latitude. Although the dust 
grains are not smoothly distributed, the approximately linear correlation between hydrogen 
scattering column density and the effective distance is expected statistically. Figure 7 
shows the correlations with the typical scale height d = 100 pc for both MRN and WDOl 
models, where Nh is only for x = 0.0 — 0.99, i.e., Nh in a; = 0.99 — 1.0 is not included 
since we believe these scattering media are localized around the sources and should not be 
correlated with the distance from the observer. The sources whose galactic latitudes are less 
than O.S'^ are also not included, because 0.1/sm(0.3°) > 16 kpc, i.e., these sources are within 
the dust layer and the real distance is less than the effective distance of d/sin{\b^^\). We 
notice that the source GX 3+1 is signihcantly below the lines. One possibility is that this 
source is actually located within the galactic dust layer and thus its true distance is much 
less than the effective distance. On the other hand if the correlations are used as distance 
indicators for sources in the dust layer, the htted value would indicate a distance of 4.8 or 
5.3 kpc (as indicated in the figure), which is in remarkable agreement to 4.5 kpc obtained by 
Kuulkers & Klis (2000) based on the Eddington luminosity X-ray bursting property of GX 
3+1. As a sanity check, we also examined the correlations with either the column density 
in a; = 0.99 — 1.0 or the total column density in a; = 0.0 — 1.0. The former correlation is 
extremely weak, consistent with the interpretation the column density in a; = 0.99 — 1.0 
is localized to each source. The latter correlation, though statistically signihcant (largely 
because the column density in a; = 0.99 — 1.0 is much smaller than in a; = 0.0 — 0.99 for most 
sources), is still much poorer than that shown in Figure 7. 


3.5. Evidence of molecular clouds around these X-ray binaries? 

In section 2, we have suggested that there should be substantial dust grains near the 
sources and these high density dust may be the molecular clouds in which the binary systems 
were formed initially. Using the same method for htting the halo prohle of LMG X-1, it is 
possible to determine the sizes of these molecular clouds. We use this method to £t the halo 
profile of Gygnus X-1 as a test and the best fit parameters are: = 0.147 x 10^^ cm“^ 

in 0.994-1.0 for WDOl model and Nh = 0.153 x 10^^ cm“^ in 0.992-1.0 for MRN model. 
Assuming the distance to Gygnus X-1 as 2.5 kpc and Gygnus X-1 is in the center of the 
assumed molecular cloud, we can calculate the radius of the cloud, twdoi = 16 pc, the 
density, pwDoi ~ 30 cm“^, and the total mass, Mwooi ~ 1-2 x IO^^Mq for the WDOl model. 
Similarly vmrn = 20 pc, pmrn 25 cm ^ and Mmrn ~ 2 x IO^Mq for the MRN model. 
These results are consistent with the typical values of molecular clouds. 

We therefore suggest that the high angular resolution X-ray halos observed with Chan- 
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dra provide evidence for the existence of molecular clouds surrounding these X-ray binaries. 
Unfortunately for most of these sources the data quality, i.e., the limited accuracy in deter¬ 
mining the halos at angles around l", does not allow us to carry out detailed studies on the 
exact sizes of these molecular clouds. Future sub-arcsecond, or even much hner angular res¬ 
olution X-ray telescopes (e.g., built with X-ray interferometry technique, Cash et al. 2000), 
will allow not only the determination of the spatial distribution of these molecular clouds, 
but also investigations of the supernova remnants of these neutron stars and possibly also 
black holes, and even the stellar winds of their companions. 


4. Conclusions 

With excellent angular resolution, good energy resolution and broad energy band, the 
Chandra ACIS is so far the best instrument for studying the X-ray halos. Although the 
direct images of bright sources obtained with ACIS usually suffer from severe pileup which 
prevents us from obtaining the halos in small angles, we can use our method to resolve the 
point source halo from the Chandra CC-mode data or grating data and get the halos in very 
small angles. 

In this paper, we analyzed 17 bright X-rays point sources using our method. The halos 
between 1.5" and 60" have been obtained. Then we use the grain models WDOl and MRN 
to £t the halo radial profiles and derive the hydrogen scattering column densities and the 
spatial distributions of dust grains for these sources. Both the WDOl and MRN model 
fittings indicate that there are substantial dust grains close to the sources. We suggest this 
as evidence for the existence of molecular clouds surrounding these X-ray binaries. 

The scattering hydrogen column densities derived from the interstellar grain model 
WDOl and MRN both have the linear correlations with the absorbing column densities 
derived from the spectral model fits (power law model, thermal bremsstrahlung or black 
body), Nh,wdoi ~ (0.720 i 0.009) x NH,abs T (0.051 i 0.013) and Nn,mrn ~ (1.156 i 
0.016) X Nn^abs + (0.062 ± 0.024), in units of 10^^ cm“^. Both interstellar grain models can 
describe the observed halo profiles adequately, and it is difficult to determine which grain 
model is preferred from our analysis. 

At present we cannot obtain dust grain spatial distributions between x = 0.0 — 0.6 for 
these sources in our Galaxy because the lack of halo intensity distribution covering a wide 
range of scattering angles. However the additional data to be collected with Chandra in the 
future for these sources studied here and especially more sources throughout the galaxy may 
allow us to probe the spatial distribution of interstellar dust media in many parts of the 



Milky Way. 
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Fig. 7.— The hydrogen scattering column density vs. distance through the galactic dust 
layer. 




